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Abstract 

The Langevin dynamics of a ci-dimensional mean spherical model 
with competing interactions along m < d directions of a hypercubic 
lattice is analysed. After a quench at high temperatures, the dynamical 
behaviour is characterized by two distinct time scales associated with 
stationary and aging regimes. The asymptotic expressions for the 
autocorrelation and response functions, in supercritical, critical, and 
subcritical cases, were calculated. Aging effects, which are known to 
be present in the ferromagnetic version of this model system, are not 
affected by the introduction of competing interactions. 

PACS numbers: 05.50. +q, 05.70.Ln, 64.60.-i, 75.10.Hk 

1 Introduction 

Non-equilibrium phenomena, as aging and violations of the fluctuation - dis- 
sipation theorem (FDT), have been attracting the attention of many investi- 
gators. A number of dynamical calculations for disordered as well as uniform 
magnetic model systems point out the occurrence of aging and violations of 
the FDT in a time evolution from a quench at high temperatures P3 12 E] ■ 
Since there are no general principles to understand and classify these dy- 
namical phenomena, it has been valuable to analyse the dynamical behavior 
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of simple, analytically tractable, model systems. In the present work, a 
detailed investigation of the Langevin dynamics of an analytically tractable 
(i-dimensional mean spherical model with competing interactions is reported. 

Spin models with competing interactions are known to display a rich 
phase diagram, with multicritical points and modulated phases. In terms of 
the temperature T and of a parameter p gauging the strength of the com- 
peting interactions, the phase diagram of the axial next-nearest-neighbor 
Ising (or ANNNIjU |3]) model displays a Lifshitz point, at the meeting of 
paramagnetic-ferromagnetic and paramagnetic-modulated critical lines, and 
an impressive sequence of modulated structures at low temperatures. The 
thermodynamic behavior of a spherical version of an Ising model with com- 
peting interactions has been originally analysed by Kalok and Obermair [§]. 
A spherical analog of the ANNNI model, with the characterization of a Lif- 
shitz point and the existence of ordered ferromagnetic and helical phases, has 
been introduced by Hornreich and coworkersjTj. The field behavior of this 
spherical analog of the ANNNI model has been investigated by Yokoi and 
coworkers 8j. With the exception of a few numerical works, as the analysis 
of an Ising model with both ferromagnetic and antiferromagnetic dipolar in- 
teractions in order to account for the behavior of ultrathin magnetic films^, 
we are not aware of analytical investigations of the dynamics of statistical 
models with competing interactions. 

In a recent paper, Godreche and Luck [2| reported a detailed analyti- 
cal treatment of the Langevin dynamics of a ^-dimensional ferromagnetic 
mean spherical model. The present work may be regarded as an extension of 
this analysis for a mean-spherical model with competing interactions. The 
particular results of Godreche and Luck are recovered. 

The layout of this paper is as follows. The spherical model with com- 
peting interactions is introduced in Section 2. The Langevin dynamics, with 
the inclusion of a time- dependent Lagrange multiplier for implementing the 
spherical constraint, is analysed in Section 3, but the mathematical details 
of this analysis are left for the Appendix. In Section 4, some comments are 
made and the main conclusions are presented. 

2 Definition of the Model 

The grand canonical partition function, 
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Sjv(/3,/x) = / exp 



-f3H({S x }) -frY, S * 



subjected to a spherical constraint, 



(1) 



-i|ln= ff(/ 3,rt = J V, 



(2) 



is the trademark of a mean spherical model. The spin variables S x G R 
are continuous, /3 is the (inverse) temperature, // is a Lagrange multiplier 
that canonically ensures the spherical constraint, and Ajy = {— L, — L + 
1, • • • ,L,L + l} d is a hypercubic lattice with iV sites. The Hamilton function 
is given by 



H({S X }) — — - £} Jx,x'S x S x /, (3) 

x,x'gAjv 



where 



x — x' = ±ej i £ {1, • • • , m} 
x — x' = ±2ej £ G {1, • • • ,ra} /.x 
x — x' = ±ej z G {m + 1, • • • , d} 
otherwise 

This exchange integral describes the whole features of the model. There are 
nearest and next-nearest neighbor interactions along m (< d) out of the d di- 
rections of the hypercubic lattice; along the remaining d — m directions, there 
are only ferromagnetic, J > 0, nearest-neighbor interactions. We assume pe- 
riodic boundary conditions along each direction (Sl+i = S-l)- Parameters 
R and S are free, but the scenario of competition takes place for S < 0. In 
the particular (ferromagnetic) case analysed by Godreche and Luck, m — 0. 
The simple spherical analog of the ANNNI model is given by m = 1 (with 
the parameter p = —S/R gauging the strength of the competition). 

The partition function can be obtained by standard proceduresjlO]. In 
the thermodynamic limit, the spherical constraint leads to the relation 



J x,x' 
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d d k 



2(27r)^-|j(ife)' 

[— 7T,7Tj a 

where f3 is written in terms of the Lagrange multiplier and 



(5) 



fi> fi c := -J(k c ) = - sup \J{k) 

^ * fcef-7r,7rl d 



where 



J{k) = 2J 



R cos ki + S cos(2fcj) + cos fcj 



i=l 



8=1 



i=m+l 



(6) 



(7) 



is the Fourier transform of the exchange integral. The critical wave vector 
k c comes from equation (JBJ). Thus, one can write 



k c = (<? c , ••• ,gc,Q, ••• ,0), 

d—m 



(8) 



where this point in the first Brillouin zone satisfies the condition (|6|) . and is 
determined by the parameters i? and 5. It is easy to see that 



' , R > and S > -R/A 



q c = < 7r , < and S > —R/4 , 



[ ±0 , S<-|i2|/4 



(9) 



where := arccos (— 

The sum rule ©, which defines the critical temperature /3(// c ), also leads 
to the lower and upper critical dimensions, d c and d, which are listed in Table 
1, in terms of d, m, and the ratio r = m/d. Note that it is convenient to 
introduce and analyse five different cases. Also, note that these ingredients 
will be sufficient for characterizing the asymptotic dynamical behavior. 
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m^O, d — m^O 


Case 1 




R 


+ 4S^O 


4 = 2 


d = 4 


Case 2 




R 


+ 4S = 


4 = 2^7 


d= ±- 

2—r 


m^O, d — m = 


Case 3 




R 


+ 4S^O 


4 = 2 


d = 4 


Case 4 




R 


+ 4S = 


4 = 4 


d = 8 


m = 0, d — m^O 


Case 5 


X 


4 = 2 


d = A 



Table 1: Lower and upper critical dimensions. 

3 The Langevin dynamics 

The dynamics is assumed to be governed by the Langevin equation, 



dS x (t) 



S 



at 



ss x (t) 



H[S x ](t) +f ,(t)J2s 2 x (t)\+Ut), 
xeA ) 

where {^xif)} is a set of random variables such that 

(&(*)) = and (Ut)&tf)) = 2T5 x , x ,5(t - t'). 



(10) 



(11) 



In contrast to the static case, the Lagrange multiplier /i is now a funtion 
of time, ensuring the spherical constraint at each time t. 

In this work, the calculations are limited to the long - time behaviour 
of autocorrelations and response functions. The analysis is restricted to the 
asymptotic expansions (for large t') of the autocorrelation, 



(12) 



fceA/> 



with t > t', where Ck(t,t') = {Sk(t)S-k(t')) is a two-time correlation in the 
Fourier space A, and the response function, 



SS x (t) 



N V ShJt') 



HO 



]y E R k(t,t'), 



(13) 



fceAjv 
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where Rk(t,t') = S(S x (t))/Shk(t'), and h is just a small perturbation. Ac- 
cording to the fluctuation-dissipation theorem, in a stationary regime these 
functions are related by the expression 



X(t,t') 



TR(t, t') 
deC{t,t') 



1. 



(14) 



If X(t,t') 7^ 1 the theorem is violated, which suggests the introduction of an 
effective temperature T/X(t,t'), larger than the heat-bath temperature T, 
and which is supposed to gauge a non-stationary behaviour of the system. 

In order to calculate the two-time functions, one may first define the 
functional 



which will be denoted by ip(t). By solving the differential equation (fTUj) and 
using the definition of Ck(t,t'), one can show that 



where (7^(0, 0) is the initial condition. For a quench from a totally disordered 
state, at an effectively infinite temperature, one should take (7^(0,0) = 1. 

The autocorrelation is obtained from the spherical constraint C(t,t) = 1 
(see (|T2|) and (J2J)), that implies, in the thermodynamics limit, and for t > 0, 




(15) 



C k (t,t') 




(16) 



o 




(17) 



o 



with 
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fit) :- 



[— 7T,7r]' 



d d k 



^2j(k)t = [j (4Ji)] 



d—m 



3 4Jt9 ^dk 



;is) 



where Iq{x) is the modified Bessel function of order zero, and 



g(k) := Rcosk + S cos(2/c) 



(19) 



corresponds to the portion of the exchange integral responsible for the com- 
petition (if S < 0). 

The convolution product in equation (JTTj) suggests a solution by Laplace 
transform, which yields 



<T-t-lOO 



£[/](*) 



-2T£[/]( S ) 



(20) 



where cr is larger than any poles of the integrand. The problem is now to 
determine if), since the autocorrelation, 



C(M') 



and the response function 



f[ t -^-)+2T fdyf (^-y )ri,j) 



(21) 



R{t,f) = f 



t-t' 



W) 



can be both written in terms of "0- 
The asymptotic behaviour of / is 



(22) 



f{t) ~ # P 



; 2J(fc c )t 



(23) 



where p labels the set p = {.R, 5, m, d}. As it is fully discussed in the Ap- 
pendix, the expressions for K p and j p , which are listed in Table 2, depend 
on the four parameters {R, S, m, d}. 
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K p 


1p 


Case 1 


2 m (8vrJ) 3 




d 
2 


Case 2 


(SvrJ)-^" 1 (H)-r(|) m ^)( Qc) -* 


2d—m 
4 


Case 3 


2 d (87Tj)"^ 




d 
2 


Case 4 


(8vrJ)^(i|)^r(f) d ^)(g c ) -* 


d 
4 


Case 5 


(8ttJ)"2 


d 
2 



Table 2: K v and 7 P . 



The behaviour of ip for large times demands the asymptotic expansion 
of £[f] for small values of e := s — 2J(k c ) > 0. With the same notation, 
r = m/d, and p = {R, S, m, d}, one can calculate the asymptotic expression 



£[/](*) 



e -ap 



Fp(-lne) 



Ax -F p (-elne) 
I Aj - F p e 



< d < d c 
d = d c 

d c < d < d , 
d = d 
d>d 



(24) 



with the coefficients given in Table 3. Note that the nth derivative 

of g, given by equation (fTTJj). which is different from g p . Also, note that 



a. 
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Case 


F P 


A 
(1 


Case 1 
a = 4=2 

V 2 


X 


a < 2 


2- (%njy l \g( 2 \q c )\-^ 


a = 2 


2 m (8vrJ)-5 |^ 2 )(g c )|-f 


2 < a < 4 


2 m (8tt J)~ 2 




a = 4 


A 2 


(7 > 4 

Li ^ T: 


Case 2 

2d-m-4 
a p 4 

g P = r (±=^) 


(SvrJ)-" (S)' r r(|) m |^)(g c )|-f 


I 4 

2— r 


(SvrJ)- 1 (f)^r(f) m |^)(g c )|-f 


7 4 

a = 

2— r 


(ftrJ)-~ (H)-r(|)- ,(% c) - 


2—r 2—r 


(svrj)- 2 (5)-r(|)^ ^)(, c )|-t 


A 8 

a = o — 

2—r 


^2 


d > — 

U> ^ 2 -r 


Case 3 

LXp 2 


2 d (87rJ)^|^ 2 )(g c )|-| 


A ^ O 


2 d (87TJ)" 1 




r\ — 9 


2 d (8ttJ)^ 


^ (2) fe)i-^ 




2 d (8ttJ)" 2 


S (2) (<7e)|- 2 




A 2 


d > 4 


Case 4 
a - d ~ A 

P 4 


(8njyH$yT(i) d \ g ^( qc )\-i 


A ^ A 

a < 4 


(SttJ)" 1 («)r(f)*bW( gc )|-^ 


a = 4 


(87rJ)-^(f)^r(f) d |^)(, c )|-! 


4 < < 8 


( 87 rJ)- 2 (frr(i) 8 |,(% c )|- 


(7 = 8 


^2 


d > 8 


Case 5 

LXp 2 

9 P = rm 


(8ttJ)~^ 


d < 2 


(87TJ)" 1 


d = 2 


(8ttJ)~2 


2 < d < 4 


(87rJr 2 


d = 4 


^2 


d > 4 



Table 3: F p , g p , and a 
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The next step is the determination of the asymptotic behaviour of i/j. The 
calculations are separated in three parts, each of them corresponding to a 
different temperature regime. 

3.1 Supercritical dynamics 

If a system is quenched from a highly disordered state (for instance, the 
system may have an effectively infinite temperature) to T > T c , the function 
ip has an asymptotic exponential behaviour 

ij>(t) ~ e t/Tp , (25) 

where r p is related to the characteristic time. This behaviour indicates the 
decay of the system to an equilibrium state in finite time. In this situation 
the autocorrelation, 

oo 

C(M') ~ C(r) =T J dyf (|) e -#, (26) 

T 

and the response function, 

R(t,t')~R( T ) = fQe-^, (27) 

depend on the time difference r only, and the fluctuation - dissipation theo- 
rem is satisfied, 

X(t,t')~l. (28) 

3.2 Critical dynamics 

In contrast to the other cases, the critical dynamical behaviour depends on 
the dimension of the system. In the following calculations, it will always be 
assumed that d is larger than the lower critical dimension d c , so that T c ^ 0; 
in other words, the occurrence of a phase transition is assumed. 
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The asymptotic behaviour of ip is given by 



' e 2J(fc c )i 



d r < d < d 



e 2J(k c )t 
hit 



d = d 



(29) 



k e 2J(fcc)t ; d>d 



which is sensitive to the dimension. 

Two time scales arise in the analysis of the critical dynamics: 

(i) In the stationary regime, 1 ~ t < f , both the autocorrelation, 



C(t,t) ~ C^T), 



with 



(30) 



C eq ,c(r) :=T e j dyf (|) 

T 

and the response function, 



-J(k c )y 



(31) 



R(t,t')~f 



-J(k c )r 



(32) 



are invariant under time translation. The fluctuation - dissipation theorem 
is satisfied with X(t, t') ~ 1. The choice r ~ 1 precludes the system to decay 
from the stationary state, which suggests the occurrence of aging for larger 
values of r. 

(ii) For 1 r ~ t' , it is convenient to define 



x := l -. (33) 



In this regime, the autocorrelation, 



11 



2K p T c 2~>P ,/1-tp a: 1 "^^-!) 1 -^ 
7p— 1 x+1 



C(t,t') ~ < 



x [(a; - 1) 



1 _i_ a£ v 
1 + T^F X 



1-7 P 



— x 



d r < d < d 



d = d 



(34) 



7 P -1 IA 



X — 1 



[X 



and the response function, 



/ 2>iO'~ 7p a;^ £ (a; - l) _7p 



d r < d < d 



R(t, ~ <( 2^K p t'~ lp (x - ir 7p \/l + , d = d 



(35) 



L 2^K p t'~ lp (x - 1 



-7p 



d > d 



show that the time translation invariance is broken, and aging effects are 
observed. 

The asymptotic behaviour of the fluctuation-dissipation ratio is calculated 
from the equations plj) and (}35|) . which yield 



X(M') 



( 7p -l)(x+l) 2 



(7 P ai+7p-2)(x+l)-2(ir-l) 

2( 7 p-l)lnt / 

2(>-D[l + (f^) 7 1ln*'-(,-l)[l-(f-i)^- 1 



1+ 



-1 \7p 



d c < d < d 
d = d 

d > d 



Note that x ~ 1 is the stationary limit, and X(t,t') ~ 1 in this case. 



(36) 



3.3 Subcritical Dynamics 

Again, the occurrence of a phase transition is required, and the calculations 
are performed for d > d c . The asymptotic behaviour of ip is given by 
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where 



(37) 



Mi := 1 - 



'eg 



T 



is the square of the static magnetisation. 

In the stationary regime, 1 ~ r <C t', the autocorrelation, 



(38) 



C(t,t')~Ml+(l-M* q )C eq AT), 
and the response function, 



(39) 



J(k c )r 



(40) 



depend on r, and the fluctuation - dissipation theorem is asymptotically 
satisfied. 

In the aging time scale, 1< t ~ the autocorrelation, 



C(M') ~ M e 



eg 



4x 



(x+ir 



2p 

2 



and the response function, 



(41) 



R(t,t') ~ K^t'^x^ (x - l)~ 7p , 
are not invariant under time translation. One may calculate 



(42) 



lim lim C(t,t') = M 2 = 1 



t— >oo V — >oo 



T 

-<■ c 



(43) 
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which is the analogous of Edwards-Anderson order parameter. This is a 
connection between the two time scales, and one can also interpolate the 
autocorrelation as 



Ax 



2 



C(t,t>)~(l-M 2 eq )C eq , c (T)+Mi 



(x + 1) 2 



(44) 



The fluctuation-dissipation ratio is 



2TK P 



( 



x+1 



) 



1+7p 



X(M') 




(45) 



x — 1 



4 Conclusions 

In this work, the Langevin dynamics of a (^-dimensional mean spherical model 
on a hypercubic lattice with nearest-neighbor (J and RJ) interactions and 
the addition of extra next-nearest-neighbor (SJ) interactions along m < d 
directions was analysed. For S < there is a scenario of competition between 
ferromagnetic and antiferromagnetic interactions, with the occurrence of an 
ordered modulated region in the phase diagram (in terms of the temperature 
T of the heat bath and the competition parameter p = —S/R). The asymp- 
totic expressions (for large values of time t') for the autocorrelation, C (t, t'), 
and the response function, R (t, t'), with t > t', were obtained, and the valid- 
ity of the fluctuation-dissipation ratio, X (t,t') = TR(t,t') /d t >C (t,t'), was 
checked. 

The addition of competing interactions does not change the qualitative 
dynamical behaviour as compared to the ferromagnetic case (case 5 in this 
work), which has been analysed in detail by Godreche and Luck[2j. The 
supercritical dynamics is trivial. The asymptotic forms of the two-time func- 
tions are translational invariant, X (t,f) ~ 1, and the system reaches equi- 
librium in a finite time. In the critical and subcritical cases, one is led to 
consider two distinct natural time scales: (i) For 1 ~ r <C t', the two-time 
functions depend on the difference r = t — t' only, and the fluctuation- 
dissipation theorem holds; (ii) If 1 <C r ~ t', in general both the autocorre- 
lation and the response function 1 depend on t and t' (instead of depending 

1 Except the response function in critical dynamics for d > d. 
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on t only). This lack of translational invariance leads to violations of the 
fluctuation-dissipation theorem, and to a system that ages with time. 
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Appendix 

The simple ferromagnetic model (case 5 in the classification of this work) is 
used in order to give a detailed account of the calculations. Results for the 
other cases can be obtained by analogous manipulations. 

Lower critical dimension 

The lower critical dimension, d c , is established by the spherical constraint 
(JSJ) at the critical value fi = /i c , 



J{k c )-J{k) = J 



d 



J 

12 



=1 i=m+l 
m d 

c 4 ^2(h- q c f + k 



i=l 



i=m+l 



Jc 3 
3 

+ • 



i=l 



£ (h - q c f + 

(46) 



where c 2 := R cos q c + 4S cos (2g c ), c 3 := R sing c + 8S'sin(2g c ), and c 4 : = 
—Rcosq c — 165'cos(2g c ). It is easy to see that c 2 > 0, and c 2 = if and only 
if R + 4S = (corresponding to q c = 0) and R — 4S = (corresponding to 
q c = tt). In these cases C3 also vanishes and therefore the fourth-order term 
becomes relevant to characterize the critical behaviour. 
For the fifth case, the (inverse) critical temperature is 
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. . 1 f d d k 1 f d d k 



J J {2*) d Y* =l k* + 0^) J (2n) d J(k c ) - J(k) 

Bs A\B S 
d ^ d 



J J (2vr) d Eti^ 2 

+ 0(5), (47) 



<5 

1 2ti% /" <i£;£; 



d ■■ - ' d-3 



JV(f)J (27T) 

where -B^ is an open ball of radius 5 centered at (0, • • • ,0), and in the last 
step the (hyper) spherical coordinates were invoked. The integral converges 
for d > 2, establishing d c = 2. In this work, h = 0(x) means that h is of 
order x or less than it; by h = o(x), it means that h is of order less than x. 

Initial conditions 

From equation (|T2"j) . the autocorrelation in Fourier space at t — t' — is 
given by 



C*(0,0) = 1 £ (5,(0)^(0))e^'), (48) 

where 



N 

a;,a;'GAjv 



<&(0)S,,(0)> = { <|<^<°>> • (49) 

for an "infinite temperature" condition. In this highly disordered situation, 
from the spherical constraint C(t, t) = 1 at t = 0, one has 



xeAjv 

= iV(^ 2 (0)), (50) 



from which (S£(0)) = 1. 

Therefore, (5^ (0)5^/(0)) = S x>x i, which yields (7^(0,0) = 1. 
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Asymptotic behaviour of / 

For large t, and choosing 5 <C 1 such that t~ 1 / 2 <C 5 <C t~ 1//4 , one shows that 
(case 5) 



/(*) 



-i d 



71 



71 



dke 



4Jt cos k 



dke 



4Jt(l-^+C(<5 4 )) 



7T 



dke 



4Jt cos fe 



D 4Ji 



7T 



3 4Jt 



O 7T 

J dke~ 2Jtk2 (1 + C(M 4 )) + i y rffce 



4 Jt cos fc 



1 v^tF 



„4J<ft 



erf(72Jt5) + C(t 3/2 5 4 )) + C(e 4J * cos<5 ) 



(8?rJt): 



(51) 



In general, one has equation (|2*3|) with K p and 7 P given in Table 2. 

Asymptotic behaviour of C[f](s) 
One should consider 



d d k 



1 



(2?r) e + 2 J(Jfc c ) - 2 J(Jfc)' 



(52) 



where e := s — 2J(k c ) > 0, and recall that fc c = and 2J(k) = 4 J Ylt=i cos ^« 
in case 5. In order to obtain the Laplace transform of / in the vicinity of 
2J(k c ) (or e ~ 0), and using the same notation as in equation ()47|). this 
expression is rewritten in the form 
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£[/](«) 



d d k 



B S 
+ 



(2tt) e + 2 J(k c )- 2J(k c )-2JJ2tikHO(5 A ) 



+ 



-7r,7r] d \_B 5 



(2tt) e + 2 J(jfc c ) - 2J(fc) 



(53) 



Changing the first term of (J53J) to (hyper) spherical coordinates, and since 
the second term is analytic in e, it is possible to write 



2JS 2 

where 



A'; 



d d k 



{-TT,TT] d \B S 



(27T) 



2 J(fc c ) - 2 J(Jfe) 



(55) 



Also, define 



(27T) 



2J(k c ) - 2J(k) 



(56) 



From equation (jSJ), it is easily seen that A\ = ^r. 

Let d = 2q. For even dimensions, and choosing e <C < 1 such that 
eln5 <C <5 2 and 5 2 <C |elne| (a possible choice is given by 5 2 = eyj In e In 5\), 
the integral in equation is 
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r> 7t£C; 1 ) ( - i), " 1 " m " m 

J m=0 v 7 

' -Ine + C(ln5 2 ) , q=l (d = 2) 

< elne + C(5 2 ) , g = 2 (d = 4) . 

k (-l) q e^lne + 0(S 2 ^) , q>2 (d > 4) 

(57) 

On the other hand, for non-even dimensions d G (0, 2), one sees that 

An analytic continuation from (0, 2) to R\Z using the functional equation 
ASH leads to 

2.J.S 2 



Therefore, 



2Ji 2 

y 
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£[/](*) 



3=1 



57T, 



+ < 



87rJ) _1 (-lne) + 0(ln5 2 ) 
8vrJn T (1 - f ) + 0(5 d - 2 ) 
87rJ)" 2 elne + C(5 2 ) 



d < 2 
d = 2 
2 < d < 4 
d = 4 
d > 4 



(60) 



One should now analyse the behaviour of Y^jLi ( — 1 ^ e '' 1 f° r e ~ 0. 
First, note that Aj (equation (JHEJ)) is finite for d > 2j. In this case 



A'- 



Aj — 



d d k 



1 



(2k) 



2J(k c ) - 2J(k) 



2vrf 



(2vr) d r(f) (2jy 



dkk d ~ l ~ 2j + 0(S 2 



A J + 0{5 d - 2] ) + 0{5 2 ). 



(61) 



For d < 2j the integral Aj diverges, and one should evaluate the asymp- 
totic behaviour of e*~ l A'j to add to (jHUjl and, therefore, characterize C[f](s) 
for s ~ 2J(A; C ) (or e ~ 0). Using cosx < 1 - ^ for x e [0, k] C R, 
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2-d 

2 



^J)-5 r (1 - §) e 



(87rJ) _1 (-lne) 



£[/](*) ~< Ai-(87rJ) _ '|r(l-f)|e 
Ai - (87rJ) _2 (-elne) 
I Ax - A 2 e 



ri-2 
2 



d < 2 
d = 2 
2 < d < 4 
d = 4 
d > 4 



Asymptotic behaviour of ^ - general comments 

In this Section the asymptotic behaviour of ijj, 



*W = 5S / 



dse s *£ 



£[/](«) 



l-2TC[f](s) 



(64) 



(65) 



where cr is larger than the real part of any pole of the integrand, will be 
evaluated. First, one should note that the Laplace transform of /, 



£[/](«) 



d d k 



(27T) 



2J(fc) 



(66) 



has a cut on [inf fee r_ 7r7r id{2J(A;)}, 2J(A; C )] in the complex s-plane. Further- 
more, C[f] is a monotonically decreasing function of s, ranging from to 
f3 c /2 (Pc/2 being infinite in the absence of phase transition). 

Asymptotic behaviour of ip - supercritical case 

Let P := {y G C : y is pole of e st C[ip](s)} and let p := sup^ gP {Re y}. By 
the monotonicity of C[f)(s), the denominator of C[ip](s) (see (|S5))) runs the 
interval [1 — |r, 1] reaching each point only one time. Therefore, in the 
supercritical dynamics, T > T c , equation (J63j) has a single pole, denoted 
henceforth by r~ x . 
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MAAA/WWWW 



Figure Al: Contour for integration - supercritical case. 



By the residue theorem, choosing the contour indicated in Figure Al, 
where ABCD is a rectangle with vertices c±iR and <j±iR such that 2J(k c ) < 
c <p < a, 



a+iR 



2ni Res e st L[ip}(s) 



dse st Cyj}{s) + e lRt / dye yt C[^}{y + iR) + 



-iR 



-R 



+ie ct J dye iyt C[^](c + iy)^ 

R 

a 

+e- lRt J dye yt C[i)}{y-iR) 



(67) 



It is easy to see that Hindoo | C[ip] (y ± iR)\ = 0. Moreover, the third 
term is 0(e ct ), which is negligible as compared with the first one (equal to 
2niif>(t) in the limit R -> oo), that is 0(e at ). Therefore, 



tp(t) ~ Res e st C[ip}(s) 



1 



4T 2 a s £[^](r- 1 



-e T p . 



(68) 
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As the temperature gets close to T c (from above), r p 1 became closer to 
2J(k c ) (hitting this point at T c ), which is one of the edge of the cut in the 
complex s-plane. Hence, in the vicinity of T c , one has 



Res e st C[i)]{s) ~ lim 

s~2J(fc c ) H 



s - 2J(k c 



e st Cm(s 



T ~ T c +, (69) 



and £[ip] can be replaced by its asymptotic formula 



emu 



Fp9v 



e- a P-2TF p g p 

F p (-\ne) 
l-2TF p (-lne) 



d < d c 
d = d r 



l-2TAi+2TF p | 9p |e Q P ' d c < d < d , 



(70) 



A 1 -F p (-elne) 
l-2TAi+2TF p (-e In e) 

A 1 -A 2 e 
k l-2TA 1 +2TA 2 e 



d = d 
d>d 



which can be calculated from (|24j) and ()65jl . If these results are inserted in 
(JBIIJ), one finds 



2J(k c ) + (2TF p g p ) 



2J(k c ) + exp 



2TF Z 



f ~ < 



2J{k c ) + 
2J(k c ) + e 
2J(k c ) + 



1 - 3L 

2TFp\g p \ V T c 



d < d r 



d = d r 



, d c < d < d 



2TA 2 



d = d 



d> d 



(71) 
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where e* is the least root of 



elne = — \— (l- ^ ) . (72) 



2TF p V T c 



The characteristic relaxation time, T eq , is related to r p by 



r" 1 =T- 1 -2J(fc c ). (73) 

Asymptotic behaviour of ip - critical case 

The simple pole of C[ip](s), which is isolated in the supercritical case, touches 
2J(k c ) (one of the edges of the cut) at T = T c . Taking the integration contour 
as in figure A2, it is easy to show that 



o+iR 

J dse st £{iP}(s) = J dse st C[^](s), (74) 

cr-iR GFEDC 

since the contribution of the paths BC and GA vanishes in the limit R — > oo. 



25 




For sufficiently large time the integral in (|74|) is dominated by the contri- 
bution due to the path FED. Therefore, substituting C[f](s) by its asymp- 
totic form (|24j) is a suitable operation to evaluate the asymptotic form of 

If d c < d < d, one has C[f](s) ~ A x - F p \g p \e ap . Hence, C[ip)(s) ~ 
-^^e~ ap and 



Fp\g P \ 



A\ 1 
F p \g p \ 2wi 



dse 



St 



s - 2J{k c 



(75) 



GFEDC 



The integral in (|73j) is a Gamma function with Hankel's contour. There- 
fore, 



^a p -l e 2J(k c )t 

? n \q„\ r(a) 



(76) 
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For d = d, C[f](s) ~ A 1 — F p {— elne), and, therefore, C[ip](s) ~ — F t \ ne . 
The integration contour is again the one shown in figure A2. Since the 
contribution due to the arc FED is zero (letting its radius to zero), one finds 



Ai e 



2J{k c )t 



2T C F P 2%i 



d( re - i7T )e tre ~ 



+ 



re l7r (lnr — in) J re m (In r + in) 
o 



d(re in )e tre 



_A ] _ 

2T C F P ( 



,2J(k c )t 



dre 



-rt 



(In 2 r + 7T 2 ) ' 



(77) 



Adopting the change of variables r — > e 7 ™*, and integrating by parts, one 



has 



A 1 


e 2J(k c )t 


2T C F P 


71 


A 1 


e 2J(k c )t 


2T c F p 


7T 



r 2 + l 



exp (-te") tan -1 r|~ + 



+nt J dr tan -1 r exp (nr — ie 7 "") | 



A 1 


e 2J(k c )t 


2T c F p 


7T 


A 1 


e 2J(k c )t 


2T C F P 


7T 


A l 


e 2J(k c )t 



2T C F P hit 



71 

2 

7T 

2 



due "tan 1 







-In ( 









due u 



7l 



7T 



In (a 



+ C(ln~ 3 t) 



du 



i In it 

1 hit 



+ C(ln" 4 t) 



(78) 



.0 

The application e~ u (l - ^f) _ function of u, is defined on the interval 

(0, oo) almost everywhere, and the integral 
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Mt) ■= I du-^ (79) 
o ln * 

is understood as the Cauchy principal value 



A{t) := v.p. / du—^. (80) 
o lnt 

In the next Section it will be demonstrated that lim^oo A(t) = 1. Using this 
result, one finds 



A e 2J(k c )t 

for d = d. 

If d>d, C[f](s) ~ Ai - A 2 e and ~ ^^TTTJ - Usin S a S ain the 

contour of figure A2, this leads to 



A 2 1 /" <= s * 

V>(*)~^r — / ds , (82) 

W A 2 2« 7 s-2J(k c ) 1 ' 

JEFGK 

with the contribution due to the path GF and DC being negligible as com- 
pared to the one due to FED. Letting its radius to zero, one sees that 



4( t ) ~ ^e 2 ^'. (83) 



The results are summarized by 
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1p(t) ~ < 



e 2J(k c )t 

F p \g\T{a p ) 

^ e 2J(fc c )i 

2T C F P lnt 

A 2 - 

±M 9J(fc c )t 



^ A 2 

Proof of lim^oo = 1 
Considering t>e, the integral 



d c < d < d 



d = d 



d> d 



(84) 



oo 

A(f) = v.p. /* du-^ = AU(t) + ^(t) + A^(t) 

i 1 lnt 



is divided in three parts such that 



(85) 



i 

lnt 

A^(t) := J du- 



i 

lnt 



In M 

lnt 



lnt / du 



lnt + I In u I 



(86) 



Therefore, 



|A«(t)| < 



lnt 



In t + In In t 



lnt 



lnt + In lnt 



due 



1 



lnt 



lnt + In lnt 



1 — e ln * 



- + 0(ln-t) 



(87) 



and lim t ^ 0O A( 1 )(t) = 0. 

The second term, A^(t), is responsible for the non zero contribution of 
the asymptotic behaviour of A(t), 
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hit lnt 



A^\t) := [ du-^ = I due~ u [1 + o(l)] = - + o(l) ; (88) 



lnt 

Tnt 



hence, Hindoo A (2) (t) = 1. 

It remains to show that the third term, A^(t), is zero for t — > oo. Sepa- 
rating in three parts, 



oo 

A< 3 >(t) := v.p. / du^—^ = A (3a \t) + A^\t) + A (3c \t), (89) 

lnt 1 " ^ 

it will be shown that each one of them vanishes in the limit t — > oo. 
The first term, 



/—U P — £l£ 

d„ I i^ = -tlnt y du^, (90) 



lnt ln * HH 



admits the bound 



1__L- l__L. 

lnt lnt 
t« 



,^, (t)j < Unt |^<_^i_y dMe -.» 



lnt lnt 
t t 



tint 



lln^fe)!' 



t 



from which lim^oo A( 3a )(t) = 0. 
The second term, 



lnt 

l- 



(91) 



du yz^ = - tlntv -P- / ( 92 ) 

In f " 



lnt 
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can be bounded by 



\A^ b \t)\ < tint v. p. / du 



- tu 



I hiu\ 



< te lntv.p. / du-—— 

J I In 



i 

" lnt 

1 
lnt 



te itt) lntv.p. y 
i_ 

lnt 
1 

lnt 

te~*( 1- ^)ln£ v.p. J 
i 

lnt 



| lnw| 



| In (1 + w) 



0?W 



= te~*( 1_ ^)mt 



and lim t ^ 00 A( 36 )(t) = 0. 
Finally, the third term, 



A (3c) (t) := 



1 



lriM 



= — tint / du 



'( 1+ Et) 

also goes to zero in the limit t — > oo, 



lnt 



ln-u ' 



i+lnT 



|^ 3c )(t)| < 



tint 



^ (1 + 



due~ tu = 



tint 



^ + 0(ln- 2 t) t 



which implies lim^oo A^ 3c ) (t) = 0, and completes the proof. 
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Asymptotic behaviour of ip - subcritical case 

In the subcritical case, where C[f](s) ~ A±—h ( s — 2J{k c )J for s ~ 2J(k c ), 

h being a continuous function for s > 2J(k c ), and /i(0) = 0. The form of h 
is given in (|24jl. It is easily seen that 



M 4 

eg 



h(s - 2J(k c 



s ~ 2J(A; C 



(96) 



where 



T 

M 2 := 1 

eq rp 



(97) 



In analogy to the manipulations in the critical case, and using the contour 
given in figure A2, one finds, for large times, that 



M 4 

eq 



dse st C[f}(s) + 0[ h(s-2J(k, 



GFEDC 



(98) 



The integral in equation ([98)1 can be written as 



1 c-^mt) - m 



M 4 

eq 



M 4 ' 

eq 



(99) 



Three auxiliary equations 

The following three equations are useful in the calculations of the asymptotic 
behaviour of the two-time functions (autocorrelation and response function): 



2T„ 



i) -L = / dte-»Mf{t). 



(100) 
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This equation is the definition of critical temperature, and it is easily obtained 
by the definition of /. 



CO 



(") ^ = j dteSfit) (T>T C ). (101) 
o 

From equation (fS5|). and r" 1 being the simple pole of C[ip](s), one has 

£if](r P 1 ) = ^. (102) 

Therefore, 

CO CO 

/ dte-*f(t) = lim / dte- st f(t) = lirn^ £[/](*) = (103) 



which is the desired result. 



in 



2T c M 2 eq 



CO 

J dte- 2i ^V(*) ( T < T c)- (104) 



o 



For T < T c , one sees that £[/](«) ~ A 1 + h(s - 2J(k c )) for s ~ 2J(k c ) + , 
where h is a continuous application of s — 2J(k c ), and therefore, h(0) = 0. 
Hence, by this expansion and from equation (fS5|). 



lim / dte~ st ip(t) = lim £ 







)2 j (fec ) + l-2T£[/]( S ) 1-2TA! 1 ; 



The desired result follows from Ai = 
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Autocorrelation, response function, and fluctuation - 
dissipation ratio 

The autocorrelation ()21jl. 



C(t,t') 



1 



i 



f(t±^] + 2T [dyf (t±l-y\^y) 



r 

f(t'+ T -)+2T J dyf (t-y + Vj^y) 



(106) 



the response function, 



W) 



(107) 



and the fluctuation-dissipation ratio, 



X(t,t') 



TR(t, t') 
d t >C{t,t 



(108) 



display different behaviours, that depend on the temperature and the chosen 
time scale. The calculations of these functions will be divided in three parts, 
each of them corresponding to different choices of temperature. In each part, 
distinct time scales, leading to the dynamical behaviour of the two-time 
function, will be considered. The notation 



will be used. 

Let if) a be the asymptotic form of ip. In other words, 



(109) 
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D >e t/r P 


T>T C 








, T = T C 


and 


d c < d < d 


n e 2j(fc c )t 


, T = T C 


and 


d = d 


Jj 3 _ e 2J(k c )t 


, T = T C 


and 


d>d 


/(*) 


, T<T C 







(110) 



Choosing an e > such that 1 <C et' <C t', one can write (jll)6|) as 



^ e J(fc c )(2t'+r) /• ft- J(fc c )(2t'+r-2y) 

+ 2T dy^——^ r ^y) + 



& + 1) 



Tp 



1 



et' 



^ e J(fe c )(2t'+r) e J " (fec)(2 ' ,+r) 
H t 1 — x 



(t' + f) 7p (t' + f) 



7p 



x / dye~ 2J(kc)v ij(y) + 2T 
o 



dyf(y)lp a [t' + g "I/ 



(111) 

Since the function = ijj(y)e~ 2J ^ kc ^ y is positive and non - increasing 

(dw(y)/dy < 0) on the real line, equation (jllljl . can be written as 



C(t,t') 



1 



et ' e J{kc)(2t'+T) 



+2T y ^/(^(V + ^-y 



(112) 
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Supercritical dynamics 

By equations (jllOj) and ()112|) . one finds 



C(t,t') ~ Q\ j——^ +2T / dyf(y)e-y^ 



(* + i) 







ef/e (J(M-^)(2t'+r)' 



-2T 



dyf{y)e-^. 



oo 

+ 2Ty dyf(y)e 



y/ T v 



113) 



(l-e)f+r/2 



In the asymptotic limit i' ~ oo, the first and third terms are negligible as 
compared with the second one, which is 0(1) (see (|l()lj) ). Therefore, one has 



C(t.t')~T I dyf (|)e 



;ii4) 



By equations (|55|) and ()107jl . the response function is 



r 



(115) 



Using equations (|1 14|) and (|115|) . one checks the fluctuation-dissipation 
theorem, 



X(t,t)~Tf 



9 m 



dyf \ 



;ii6) 



36 



Critical dynamics 

In the stationary regime (1 ~ r C t'), one may rewrite (jll2j) as 



C(M') 



et'+r/2 

f J(k c )(2t'+r) \ r 

O I -4 ^— \+2T r I dyf{y)x 



(*' + !)' 

(l-e)t'+r/2 



et'+r/2 



(117) 



which may be convenient for calculating the analytic asymptotic form of the 
autocorrelation in this regime. 

Critical dynamics (d c < d < d) 

By (jllOj) and ()117j) . in the stationary regime, one has 



C(t,t') ~ [(t' + r)t'Y 



- - p 



o 



et' 



+ 



2T r 



et'+r/2 



(t' + r/2)>7 + 



dyx 



(l-e)t'+r/2 



:f(y)e- 2j ^ + 2T c K p 



z/> (f + 1 - y) 



2-7p 



Performing the change of variable y — > 1 /y in the last term, one finds 



et' 



+ ~ : 9--V C eq c \j) 



;ii8) 



where C e q iC is defined as 
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G,„(-) / / dyf(^)e- J ^y 



(119) 



Restricting the range of e > to be such that 1 < et' < (f)min{i,2 7p -2} 
leads to 



C(t,ti) ~ C^ir) . 
Since the response function is 



(120) 



2 / y e 2j(fe c )t /fi-op 



/(- 



27 (i + f)^ 



/ (2) e- J(fcc)T , 



(121) 



it is also a function of r only, and the fluctuation - dissipation theorem, 
X(t,t') ~ 1, is asymptotically obeyed. 

In the aging regime, from ()112j) . one has 



C(t,t') ~ [(t , + r)tfp St 







et' 



(l-e)t'+r/2 



X 



2T c K p x 



[(t' + r)t']^ O 



2T C K P 



y-yp (V + \ - y) 
et' 



2-7p 



( 7p - 1) (tf + I) V T 



(* + I) 



7p 



(122) 
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Restricting e > to be such that 1 < e(t' + r/2) < (t'/r)^- 1 ^' + r/2)>, 
one finds 



CM) .g^!M. (123) 

7p — 1 x + 1 

From equations (f2"3|) and (jl()7[) . one calculates the response function, 



,2J(k c )t' /f' 1 ~ ap 
e 2J(k c )t lf}-a p 



f (I) e'W? 

2^K p t'~ Jp (x - l) -7 " x^, (124) 



1— a B 



and the fluctuation-dissipation ratio, 

Xtt, ~ 7 (7p-1)(^ + 1) 2 (125) 

1 j ( 7p x + 7p -2)(x + 1)-2(x-1)' 1 ; 

which comes from equations (jl23J) and ()124jl . 
Critical dynamics (d = d) 

In the stationary regime, 1 ~ r <C t', for d — d, one should proceed analo- 
gously as was done in the stationary regime of case d c < d < d. Therefore, 
by choosing e > such that 1 < d! < (i! + r/2) 7 f/ln(t' + r/2), one finds, 
by equations and ()117j) . that 
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C(t,f) ~ y/ln(t' + T)\nt 

et'+r/2 



o 



et' 



2T C 

+ ; rX 



(*' + i) 7 V 111 (*' + !) 

(l-e)t'+r/2 

^/(2/)e' 2j(fcc)y + 2T c iT p 



dy 



etf+r/2 



y> In (t' + | - y) 



v/ln(t' + r) lnt' 
/ 







2T C 

+ : ; —X 



X 



\ 



dy/ (y)e~ 2 ^ - 



K„dy 



yip 

/ (l-e)+r/2 \ 

}_ [ dy_ 

et' J y~iv 

\ ef+r/2 / 



+ 



~ C e g )C (r). 



Furthermore, by equation (jl(J7|) . it can be shown that 



(126) 



R(t,t') ~ /( 



r\ / e 2j ( fcc )*'/lnf 



2/ V e 2J (^V In (f + r) 

T\ e -H k <:)T 

2 



1 + 



ln(l+r/tQ 
lnt' 



T 



ensuring that X(t,t') ~ 1 in the stationary time scale. 

In the aging regime, 1 < r ~ f, from (jllOj) and (jll2j) . one finds 



(127) 
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C(t,0 ~ Vint In f 



= Vint In t' 



Vint In t 



O 



o 



o 



6? 



(t'+f)> 

et' 

(t^TiT 

et' 



(l-e)t'+r/2 



2T c K p 



dy 



y> In (f ' + \ - y) 
dy 



et' 



(t' + r- y y p lny 



x 



/du 
(t' + l-t'uVUl + ^ 



2T c K p t' 
H ^—r-x 



hit' 



(128) 



where the change of variable y = t'u was performed in the last step. Since 
the condition et' 3> 1 is equivalent to 1 3> — In e/ lnt' = | lne/ lnt'|, one finds 



C(t.t') ~ Vint lnt' 



O 



et' 



(f + i) 



7p 



+ 



2T c K p t' 



du 



lnt' J (tf + z-tfu) 



(129) 



which leads to 



C(t, t') ~ ^^(t')^> [(* - I) 1 "* - (s + l) 1 "*] y 1 + • (130) 



7 P 

The calculation of the response function is simpler, 



fc-,2W' I (131) 
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Therefore, the fluctuation-dissipation ratio is 



X(M<)~ 2fa-l)l. lt ' 

a^-ijii+dsriinf-tx-ijfi-cgir- 1 ] 

Critical dynamics (d > d) 

For d > d, by equations (jllOj) and pi2j) . the autocorrelation is 

; v (l-e)t'+r/2 

c(t,t') ~ of ^ j+2r c / d2//(y)e ~ 2j(fcc)y 

2 

C eg , c (r) , (133) 



since C eqfi {r) is 0(1), which is much larger than other terms in the asymptotic 
limit. 

The response function is obtained from equations (jSljl and ()107J) : 



R(t,t') ~ f (134) 



The asymptotic expansions of the two-time functions in the stationary 
regime (for r ~ 1) are given by equations (|133|) and (|134|) . In this case, one 
also finds that X(t,t') ~ 1. On the other hand, in the aging scenario, for 
t>1, using (|112jh these functions have the following asymptotic behaviour: 



(l-e)f+r/2 



(t' + ^J c p J y* 

2 

-^^(t') 1 " > Ux - 1) 1_> - (z + 1) 1_> ] (135) 
1p - 1 



and 
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R{t,t') ~ 2>K p r-> = 2>K/ P (x - 1) 



-7p 



(136) 



In this situation, the flutuation-dissipation ratio is violated with 



X(M') 



1 + (^) 7P ' 



(137) 



Subcritical dynamics 

As in the critical dynamics, the two characteristic time scales (stationary and 
aging) are also present. 

In the stationary case, 1 ~ t < t', from equations (J2*3|) . (jllOj) and (jlll|) . 
one has 



C(t,t') ~ M e 4 ? [(f + r) f 

(l-e)t'+r/2 

2T 
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where the equation (jl(J4j) was invoked. By choosing e > such that 1 <C 
(t') ^p- 1 < et' < t', one has 



C(M')~M;+(1-M;)Qt). 
Using equation the response function, described by 



(139) 



in both time scales, behaves as 



(140) 



7p 

J2(f, - / (0 (l + £) 2 = / (0 e-*^ + (141) 

in the stationary case. It is not difficult to see that in this case the fluctuation- 
dissipation theorem is valid, with X(t,t') ~ 1. 

In the aging regime, 1 < r ~ f, from equations (|2*3Jl and (}9T?j) in ()106j) . 
the autocorrelation can be written as 



2 
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Taking e > such that (f + r/2)r( 1 ">)/> < et' < one has 



/ 



(142) 



C(t,t') ~ M ( 



eq 



Ax 



.(x + 1) 
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2ll 



(143) 



The calculation of the response function is simpler: 



T/~ J(k )t 

R(t, t') ~ pC I e- J ^ T x^ = K^n'^x^ (x - l)~ lp . (144) 
I2J 



From equations ((T43j) and (|144|) . the fluctuation-dissipation ratio is given 
by the asymptotic expression 



X(M) ^Mf! ■ (145) 

Technical note 

From equations (jlOfij) and ()17|) . with t —> t' + t/2, one finds another form for 
the autocorrelation, 



2T J dyj(t' + T --y)i,(y) 



2 

ij (V + L^-2Tj dyf(y)ij (t' + ^-V 



(146) 



It is possible to use this expression for calculating the asymptotic forms 
of autocorrelation. One then recovers the same asymptotic results that have 
already been reported, with the exception a discrepancy in the critical dy- 
namics for d > d, which will not display any aging behaviour. Equation 
()146|h however, involves strongly varying terms, which may even change sign, 
and whose asymptotic behavior may turn out to be much more difficult to 
analyse. 
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